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The quark-gluon plasma (QGP) can be explored in relativistic heavy ion collisions by the jet quenching sig- 
nature, i.e. by the energy loss of a high energy quark or gluon traversing the plasma. We introduce a novel 
QCD evolution formalism in the leading logarithm approximation, where normal parton radiation is interleaved 
with scattering on the plasma gluons occuring at a similar time scale. The idea is elaborated in two approaches. 
One extends the DGLAP evolution equations for fragmentation functions to include scatterings in the medium, 
QQ ' which facilitates numerical solutions for comparison with data and provides a basis for a Monte Carlo imple- 

' mentation. The other approach is more general by including also the transverse momentum dependence of the 

, jet evolution, which allows a separation of the scales also for the scattering term and provides a basis for analyt- 

04 ' ical investigations. The two approaches are shown to be related and give the same characteristic softening of the 

< I ' jet depending on the temperature of the plasma. A substantial effect is found at the RHIC energy and is further 

PL(, enhanced at LHC. Systematic studies of data on the energy loss could, therefore, demonstrate the existence of 

■ the QGP and probe its properties. 

^ : 

(N . 

■ I. INTRODUCTION 

' , The quark-gluon plasma (QGP) is expected to be formed in high energy collisions of nuclei when the energy density is large 
enough. As the plasma expands it cools down and undergoes a phase transition to ordinary hadrons. Different observables have 
(— I been suggested to provide a signal for the existence of the plasma. One of the most promising signals, that can also reveal 
detailed information on the properties of the QGP, appears at present to be jet quenching, i.e. the energy loss of an energetic 
quark or gluon traversing the plasma resulting in a jet of hadrons, or single hadrons, with reduced energy IH 0> Htj- Thus, one 
may consider a normal hard QCD parton-parton scattering event to take place in the environment of the QGP which has been 
produced in the nucleus-nucleus collision. The QGP acts as a background with which the hard-scattered partons may interact. 
QQ Such jet quenching phenomena have indeed been observed at RHIC (^.'^JdJJJS], but it remains to be demonstrated that they 

CSj are really caused by a QGP. This requires a better understanding of the QCD processes involved and this paper contributes to 
CO this. Once a proper framework has been established, one can use the modified jet properties as a trustworthy signal to claim 
discovery of the QGP, and also as a tool to study its properties. 

Until recently, theoretical efforts have concentrated on medium-induced gluon bremsstrahlung as the dominant energy loss 
00 mechanism jSrafnUil 111 111. However, collisional energy loss caused by multiple scatterings in the medium has also been 
considered lllSl Ig, 1711 and demonstrated fi6l\ to contribute significantly to the jet quenching effect observed at RHIC. 
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Although these scatterings occur over the few-fermi extended plasma, they overlap in time with normal perturbative QCD 
(pQCD) gluon radiation (the parton cascade) which is time dilated in the QGP rest frame. Therefore, such radiation and interac- 
tions with the plasma must be treated in a common framework. This is the objective of this paper 
^ . The energy available at RHIC limits the transverse momentum of the hard-scattered parton such that one mainly observes 
■ ■ ' single high p± hadrons, rather than complete jets. With heavy ion collisions in LHC, however, the energy is substantially larger 
and should give access to the study of jets and how their properties are affected by the underlying interaction with the plasma. 
The larger momentum transfer should also establish a safer ground for a perturbative QCD treatment. 

In Section ini we introduce scattering with the plasma into the familiar jet evolution DGLAP equation ifTsl [T9l l2^. In this 
way the scattering can be solved numerically without further approximations and this evolution equation provides a basis for 
development of a modified parton shower Monte Carlo algorithm ll22ll . In Section |III1 we start from a more general evolution 
equation taking also the transverse momentum into account and introduce the scattering. We show how this modifies the evolu- 
tion for the dominant small momentum transfer scatterings. Using appropriate approximations, we obtain a medium-modified 
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evolution equation which is suitable for analytical investigations. We also show that by integrating out the transverse momentum 
dependence the medium-modified DGLAP equations of Section is recovered. In Section |IV| both approaches are used to 
obtain numerical results on how the fragmentation functions are affected by this jet quenching mechanism. SectionlVjconcludes 
with a discussion of our results. 



II. BASICS OF JET EVOLUTION IN THE QUARK-GLUON PLASMA 

A quark or gluon that emerges from a hard scattering can be virtual {p^ ^ m}) and therefore it will radiate successively to 
reduce its virtuality and become on mass-shell. This leads to a parton shower and scaling violations in the jet fragmentation 
functions as described by the DGLAP equations 1 18, 19, 20] and other more sophisticated evolution equations [21]. 

The formation of this parton shower does not happen instantaneously, but needs a certain time. In the lab frame, i.e. in the 
rest frame of the nucleus-nucleus collision and the plasma, the lifetime of a virtual parton can be estimated from the uncertainty 
principle XoheE / Q^, where E is the parton's energy and Q its virtual mass. In the case of a parton shower, however, the relevant 
time is the time a virtual state needs to evolve in virtuality from Qio Q + dQ, 



E dQ^ 



The time a parton needs to reduce its virtuality from the starting scale ~ E ~ 100 GeV to the hadronisation scale Qq ~ I 

GeV can then be estimated asT = ft/Tw-^ — 4. Thus, even though a high energy parton with large virtuaUty will reduce its 

Co 

virtuality rapidly, the overall lifetime can be considerable and of the order £■/(! GeV) ^ or several fermi. The formation time T; 
of the quark-gluon plasma is estimated to be of the order 0.2 fm at Rmc and likely less than 0.1 fm at Lhc, but the plasma 
lifetime Tpi^^ma — '^i{Ti/Tc)^ may be long due to the cooling from the initial temperature to the critical temperature. Taking 
into account these two time scales of the parton and the surrounding medium we arrive at the conclusion that the parton shower 
evolution overlaps in time with the plasma phase. Therefore, energy loss through radiation and scattering has to be treated 
together in a common framework. 

In this paper we present a study of a parton cascade that includes scattering off the partons in the QGP as a modification to 
the well known QCD evolution in vacuum. This is achieved by adding a scattering term in the DGLAP equations. Modified 
DGLAP evolution equations due to medium-induced radiation has been considered in ifTsIl for the QGP and for deep inelastic 
scattering off nuclei in |23]. The literature on jet quenching in the QGP also deals with effects like energy loss of on-shell and 
off-shell partons and modified hadronisation at Qq |9, 10, 11]. Our approach differs from these by explicitly handling off-shell 
partons and the evolution of the shower in the presence of the medium. 

We assume here that the hadronisation at the infrared cut-off scale Qo is unchanged and use standard vacuum fragmentation 
functions, since the plasma is essentially gone by the time the parton has evolved down to this low virtuality. Furthermore, 
we are here primarily interested in the high-px particles that stand out from the underlying background of hadrons from the 
hadronisation of the plasma and such higher energy hadrons have a longer formation time in the lab frame due to their larger 
Lorentz 7 factor. In this and the next section we concentrate on developing the formalism for the evolution at the parton level 
and thereafter add on hadronisation through a well established parametrisation. 

Let us consider a hard process where a parton with high virtuality and high energy is produced. Due to parton splitting and 
parton scattering a parton shower emerges from this source. We define the fragmentation function Dj{x,Q^) as the probability 
density for the hard parton / with virtuality to fragment into a parton or hadron j which takes a fraction x of the initial parton 
energy. In this first study we restrict our parton level treatment to a purely gluonic system, which should give the dominating 
behaviour and most essential results. The indices on the fragmentation functions can then be dropped and the formalism becomes 
simplified also in other aspects. 

In the shower evolution, an arbitrary intermediate gluon with a fraction y of the original energy E loses a fraction z through 
radiation or scattering, and becomes an outgoing gluon with energy z{yE). The parton cascade occurs through the normal 
splitting process g ^ g + g with differential probability 

splitting {Z,Q )= Pi^)^^-^ (2) 



which is of first order in and written in terms of the unregularised splitting function 

z 1 — z 

The scattering process g + g^g + gis similarly given by the differential probability 



■+z(l-z) 



(3) 



d^scat {z,yE,Q^) = ^n,{T)dzdT = aUQ^)K {z,yE,Q^) dz^ (4) 
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as the product of the scattering cross-section, the density of target gluons in the plasma and the time (or equivalently distance) 
the gluon traverses the plasma. We introduce the scattering function K in analogy with the splitting function to combine gluon 
splittings and gluon scatterings in a generalized DGLAP equation for the QCD evolution of a jet in the plasma medium 



X 

'Az{^^P{m^,Q') + a^{Q')K{z,xE,Q^)D{x,Q') \ (5) 



where in the first line ^ —y, i.e. the intermediate energy fraction y is integrated out via / dy . . .5{y — x/z). The negative 
contribution in the second line is present when using the unregularized splitting function instead of the regularized one with the 
'plus prescription' ||24|] . This formulation gives a clear physical interpretation, namely that the first line in Eq. ^ gives the gain 
of a gluon at x which was at yE before the splitting or scattering, and the second line gives the loss of a gluon with energy xE 
before the splitting or scattering. We also note that the two loss terms can be thought of as ensuring that probability is conserved. 
For the splitting loss term it leads after integration of the DGLAP equation to the so-called Sudakov factor, whereas from the 
scattering loss term one gets a factor ^ exp (— / ^tigdzdz) corresponding to the attenuation of a beam. 

To obtain the scattering function K, we develop the different parts of Eq. (|4]l. For the g + g ^ g+ g scattering cross section 
we use the leading order pQCD formula 

da 9/^9,9 2n 

which has been regularized by the Debye mass /Xd of the plasma gluons which screens the interaction such that this expression 
can be applied down to |f | = 0. Furthermore, the scale in a, is chosen as the gluon virtuality and is thereby in the perturbative 
region above Qq ^ 1 GeV^. This is as for the usual DGLAP parton splitting and describes the weaker interaction of a more 
virtual parton having a smaller gluon cloud. We use the one loop approximation 

for the running strong coupling with A = 0.25 GeV. The cross section in Eq. (|6]l implicitly assumes that the virtual mass of the 
fast parton does not change in the scattering, and applies for a fixed centre-of-mass energy, i.e. fixed yE. For a scattering centre 
at rest with mass ms, the squared momentum transfer is then 

t = 2m,yE{z-l) (8) 

The scattering is on gluons in the relativistic plasma of temperature T with density 

«.(n = ^C(3)r3 (9) 

where ^ is the Riemann Zeta function with ^(3) ~ L20. We neglect the motion of the scattering centers, which is a good 
approximation at large jet energies, i.e. xE ^ T. For small x this is not necessarily the case, but in this region the DGLAP 
equations should also be improved by soft gluon resummation for an accurate description of the shower evolution. The Debye 
mass °^ T \25] and the mass of the scattering centre is taken as = pLu/^/2. 

Including also dx from Eq. (dJ we obtain the scattering function as a function of the gluon energy S', 

where S — yE ^ ^E in the gain term of Eq. (|5]l and S" = xE in the loss term. 

This K contains the essential dynamics of the scattering, in analogy with the splitting function P in Eq. (jS). For large 
momentum transfers, i.e. — f — 2m^S{l — z) 3> jUq, the scattering function is suppressed as A' 1/(2^(1 — z)^) as a higher 
twist process. For small momentum transfers, on the other hand, the scattering function K ~ /Q^- Thus, K may compensate 
for the additional factor Us in the scattering term as compared to the splitting term in the evolution Eq. Q for large S and/or 
large x. 



X 9 \ ,X 
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The jet quenching effect is then given by the difference between the scattering gain and loss terms in Eq. Q, i.e. 



dz 



S{x,Q') - l-a^{Q^)K\z;-E,Q']D{-,Q')~ I dzaf{Q')K{z,xE,Q')D{x,Q') 



(11) 



For leading particles in the jet, i.e. for x close to 1, and for large energies, x{l — x)E ^ ms, the gain term in Eq. ( fTTT i can 
be Taylor expanded (with l/z= 1 + (1 - z) + . . . such that D{x/z, Q^) = D{x, Q^)+x{l-z) ^"''^f'^ + . . . and jiT (z, f £, Q^) = 
K{z,xE,Q^) + ...) giving 



S{x,Q^) = JdzaHQ^)K{z,xE,Q^){l-z)(D{x,Q^) 



,2^ ,_ jDix,Q^) 



dx 



dza^{Q^)K{z,xE,Q')D{x,Q'). 



(12) 



Performing the integral over z then gives 



S{x^l,Q') 



4ms ^2 



D{x,Q^)+x- 



dD{x,Q^)\ m^+x{l-x)E 



dx 



In 



(13) 



where f{x) is independent of the energy for x{l — x)E ^ m^. Thus, for large energies and large x such that, ntf^/E ^ I— x<^ 1, 
the effects of the scattering will grow logarithmically with the energy as a result of the integration limit z — x corresponding to 
the maximum momentum transfer squared |f | = 2m^Ex{ \ —x) in the scattering gain term. 



III. TRANSVERSE BROADENING IN JET EVOLUTION 



The parton transverse momenta can be explicitly taken into account in the jet evolution. Such an evolution equation in the 
vacuum was first derived in leading logarithmic approximation in ref. ll2Tll . Recently two papers ll26ll27ll have appeared which 
apply the same framework of transverse-momentum-dependent QCD evolution equations to deep inelastic scattering. Here we 
generalize this to the case of jet evolution in the quark-gluon plasma by introducing scatterings as in the previous section. 

Transverse momentum of the tagged parton is an important tool to diagnose the history of the parton in the medium. In deep 
inelastic scattering in nuclei one can trace parton scattering and from the mean p\ read off the length of the trajectory before 
hadron formation. In heavy ion scattering the jet cone related to the transverse momenta of partons will also reflect the scattering 
history. 

As before we restrict ourselves to gluon splitting and gluon scattering, since gluon dynamics dominates these processes. The 
"fragmentation function" Q^,p\_) gives the probability for an initial gluon to convert into a gluon with momentum fraction 
X, virtuality and transverse momentum /7x in the course of splittings and scatterings which give the gluon kicks. The 
evolution equation for this multiple differential function has the following form, when we generahze it to include scattering in 
the medium 

dD{x,Q\p^) ^ 

^jyjpiz,aM'))^s(z{l~z)Q'-f-qiyi^j,Q\p^--^^^^ 

(14) 

Both gains and losses from the evolution are here included by the use of the regularized splitting function P{z) (instead of 
the unregularized P{z) in Eq. (|5])). The first, integral term is the normal evolution in the vacuum. This accounts for the mass 
constraint z{i— z)Q^ = Qq/ 4- + q^cirising in the splitting with momentum fractions z and 1 — z. The transverse momenta appear 
together with longitudinal momentum fractions to guarantee boost invariance. After integration over transverse momentum one 
obtains the standard DGLAP equation. 

The second term S{x,Q^,p±^) accounts for the scattering on gluons in the plasma with temperature T and the gluon density 
ng{T) in Eq. (|9]). The scatterings change the transverse momentum of the leading fast parton by giving q± kicks, but are here 
assumed to not change the mass scale or virtuality of the fast parton. This is strictly only true for small momentum transfers, i.e. 

smaU angle scattering. The time scale for interaction with the plasma dT = (Eq lD depends on the virtuahty. As before, 
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this gives the connection between the Q^-evolution time scale and the scattering time scale and shows that the parton cascade is 
interleaved with scattering in the plasma. Combining these ingredients we obtain 



2 E /■' /■ 2 

S{x,Q ,P_l) = -^ng{T) J dy J d qj_ 



da 

X 



x{yD{y,Q\p^-yq^)~xD{x,Q\p^))5(^y~x~^^ (15) 

with a gain term for scattering into the considered p± bin and a loss term for scattering out of it. Both are weighted with the 
parton-parton differential cross section. The scattering partners of the fast parton are "quasifree" gluons with mass ms in the 
plasma. Scattering makes the fast parton lose energy which is absorbed as recoil energy ^^/(2ms) by the plasma parton. In 
contrast to the previous section we here restrict to allow only soft scatterings, i.e. g'^ ~ /i^, such that Ax = y — x = q^/ [2m^E) 
is small. 

In this paper we restrict ourselves to the px-integrated parton fragmentation function D{x,^) = J d^pj^D{x,Q^ ,pjj. The 
resulting -integrated equation for splitting and scattering is then 

where 

results from a Taylor expansion of Eq. ( fTSl ) in the small parameter Ax. This equation takes into account both the gain term 
and the loss term in a form which is suitable for analytical investigations to understand the effects of the scattering term. This 
evolution equation can be solved by the standard Mellin transform |28], d{J, Q^) = Jq dzz^^^ D{z, Q^), as will be investigated 
in a forthcoming paper. 

Here, we note that this scattering term includes the parameter q = nga(q\), which is well known from other Eikonal ap- 
proaches to energy loss |29]. The gluon-gluon cross section of Eq. ^ is of order a ~ 9naj{Q^)/{2ji^) and the mean transverse 
momentum squared approximately equals the Debye mass squared {q\) ~ ji^ where the Debye mass is determined by a self- 
consistency equation [25]. The prefactor in the scattering term Eq. ( fT6l ) can then be written as 

nJT)a(q]) T,T^ 36C(3) 

'\ ' ^ea^{Q^)-^ where £ = ^^^2— (17) 

2ms g Q n IJ-D 

For these soft scatterings one therefore finds S <^ /Q^ (recalling that /io <:><: T), which show that large (f' suppresses the effects 
of collisions compared to gluon splitting. This higher twist effect is, however, compensated by a large derivative in Eq. ( fTSI l due 
to the steeply falling fragmentation function D at x — > 1 . 

These results are similar to those obtained from the scattering function K in Eq. (fTol i when used in the gain and loss terms 
of Eq. |5] This is no accident since the general evolution equation ( fT4l l can be used to derive the modified DGLAP equation for 
evolution in the medium that was obtained in the Section |ll] To demonstrate this one starts from the scattering term in Eq. (fTSl l 
and first averages over the direction of ^x assuming that has no dependence on the azimuthal angle. Next, by integrating 

CI C/j^ 

out the px^ dependence one obtains 

where the differential cross section is different in the gain and loss terms. By performing the change of variables from ^x to z, 
with the identification —t = q\ = 2mf,Ey{l — z), and using the explicit form of the cross section one obtains for the gain term 

Scain = ^Ina^iQ^) C dy ( dq\ ^ , ^ ,^. yD{y,Q^)5 (y- - 



Q2 2'-^K-^)J^ ^±(^2+^2)2^ 2ms£ 

En 9 [ 1 

-/-7ra2(g2) / / dz2m,yE- — , , ,^. yD{y,Q^)5{yz-x) 



2 ''^'J, 'J {2m,yE{\-z)+lllY' 
9 2/^2n dz^ f X 1 



^_-.atmi -2.n. y-E^ ^_^_^_^z,| ) ^ (19, 
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Using Eqs. (|4|i and ( fTOb finally gives 



z 

for the gain term. The loss term is treated in the same way, giving 



Scain {Q')K (z,—,Q']D['-,Q'] (20) 



= J dzafiQ^)K{z,xE,Q^)D{x,Q^). (21) 

Subtracting the loss term Eq. ( 1211 1 from the gain term Eq. ( |20l l. and including the -integrated terms from the splitting in Eq. 
(fl4] i. one obtains the medium-modified DGLAP equation (|5]l. 



IV. NUMERICAL RESULTS FROM THE EVOLUTION EQUATIONS 

The formalism developed above applies for a gluon traversing a plasma of density ng{T) depending on temperature T. To 
solve the evolution equations for a plasma which expands and cools down with time is quite complicated. A realistic treatment 
of such an expanding plasma in which the traversing parton radiates and scatters requires a Monte Carlo simulation approach, 
which is presently being developed [22]. In order to obtain some numerical results to illustrate the main effects to be expected, 
we here simphfy by considering a plasma with constant temperature T — 500 MeV and corresponding density given by Eq. (|9|l, 
representing the average properties of the plasma. At this temperature /Xd ~ 1-97 GeV, obtained from fio °^ T in li25H . This will 
not give results that can be applied to describe heavy ion collision data, but should illustrate the main qualitative effects. 

The theoretical formalism above has focussed on the parton level processes. In order to compare with observable hadrons, 
one has to account for hadronization. As justified above, one may use a standard parametrisation for the fragmentation function 
at Go ^ 1 GeV^ obtained from data in the vacuum case. We use the KKP parametrisation [30.1 for the hadronisation of gluons 
into charged pions 

Dl{x,Ql) - 3.73x-<'-^42(j _^)2.33 (22) 

at the scale Qq = 2GeV^ in our evolution equations. 

The resulting gluon-to-pion fragmentation functions are shown in Fig. [T]for both the normal vacuum case and for the case 
of a plasma using our medium-modified DGLAP equation. The overall behaviour is similar. Already in the vacuum case, the 
evolution caused by a higher depletes the high x region and enhances the low x part, since the parton cascade means that 
energy is shared among more and more partons. The scattering term acts in a similar manner and "transports" partons to a lower 
momentum fraction. Our scattering formalisms account for the energy loss of the gluon traversing the medium, but not the fact 
that the struck scattering centers obtain energy that may contribute to the jet energy. Such struck partons are expected to interact 
further in the plasma and contribute to the underlying activity in the event, including possible minijets, but they will not produce 
leading hadrons in the high-px jets which is the focus of this study. 

A more detailed comparison of the fragmentation function in the medium and in the vacuum is given in Fig. |2] by their 
ratio, which shows a moderate suppression of leading particles in a 10 GeV jet and a strong one in the case of a 100 GeV jet. 
The stronger quenching effect in the latter case depends not only on the longer evolution path in In Q^, but also on the energy 
dependence of the scattering, as illustrated in Eq. (T3[ . Thus, collisional energy loss is certainly non-negligible, even using a 
perturbative cross section, and can become quite large, as also indicated in an earlier study [16J based on a phenomenological 
model for soft colour interactions. The total energy loss will be even larger when additional medium-induced radiation is taken 
into account [^.mill. 

This ratio of fragmentation functions in medium and vacuum corresponds to the nuclear modification factor Raa at a fixed 
jet energy. The Raa measured in heavy ion collisions at RHIC |8)] is for leading hadrons with p± up to ^ 20 GeV, and thereby 
integrates over jet (parton) energy. For a direct comparison with RHIC data, one therefore needs to fold this kind of fragmentation 
functions for both gluons and quarks with the initial momentum distribution of gluons and quarks from the hard scattering. A 
detailed comparison to data also requires a proper treatment of the finite-size expanding plasma. Taking this into account is 
beyond the scope of this paper, which does not aim at giving a precise prediction to be compared with experimental data but 
rather to develop a formalism and show its basic physical effects. Nevertheless, our approximation with a stationary plasma in 
the numerical solution of the evolution equations is reasonable to illustrate the magnitude of jet quenching effect. To see this 
one has to consider the time needed for a parton shower to evolve down to Qq. This depends, however, on the topology of the 
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Figure 1 : Tiie fragmentation function of gluons into charged pions obtained when solving Eq. (O for different initial virtualities and with 
the boundary condition at the perturbative QCD cut-off 2q = 2 GeV^ given by the KKP parameterization in Eq. | |22I ( for hadronisation. Left 
without scatterings, i.e. vacuum case with conventional DGLAP evolution. Right with medium-modified evolution including scatterings of 
gluons, with initial energy E = Qmax, in a gluon plasma of fixed temperature T = 500MeV compared to the vacuum case. 




Figure 2: Ratios of fragmentation functions in medium and vacuum, D^„gj{x, Q^) / ^,g^{x, 2^), of gluons into charged pions for initial gluon 
energies E = Qmax = 10 and 100 GeV (cf. Fig. [TJ. The medium fragmentation function is calculated using the modified DGLAP evolution, 
Eq. HJ, including scatterings in a gluonic plasma of T = 500 MeV. 



shower. In the evolution equation all possible radiation patterns are weighted with their probabilities and summed, making it 
impossible to keep precise track of the time during the evolution. Based on Eq. the mean time for evolution can be estimated 
s& E /Qq—1 /£, which means that a 100 GeV jet developes over lOfm to reach = \/2GeV. A finite size medium can then 
be included by defining two cut-off scales, Q\ at which scatterings stop and 2o at which splittings stop. The suppression factor 
is then found to vary roughly linearly with the path length L~ E /Q\— \/E in the plasma. For instance, a medium of length 
L = 7fm (corresponding to 2j « 3 GeV^) instead of L = °a results in a suppression factor of 0.57 instead of 0.38 for leading 
particles six = 0.8 in a 100 GeV jet. The jc-averaged suppression factor is 0.75 fori = 7fm instead of 0.64 for L — oo, ie. overall 
only an approximately 15% difference. We note that these differences are of similar magnitude as the difference between our 
two forms of the scattering term shown in Fig.[3]below. 

Furthermore, the recoiling scattering centre would also have to be taken into account. These complications can best be 
studied using the Monte Carlo framework under development [22]. Already at this stage one can, however, note that there is 
an overall general qualitative agreement, but the effect in data is much stronger than the suppression of a 10 GeV jet shown in 
Fig.|2] This indicates the need for adding medium-induced radiation and raises the question whether the medium affects also the 
non-perturbative hadronisation. 

The strongly enhanced effect of collisional energy loss at higher jet energies provides interesting prospects for heavy ion 
collisions at LHC, where individual hadrons and jets with up to ~ 100 GeV will become available. The curve for Q = 100 
GeV in Fig. |2] should not be taken as a precise prediction, in view of the uncertainties mentioned. It does, however, indicate a 
substantially increased effect of energy loss due to collisions. This is due to the enlarged phase space of momentum transfer 
squared |f |, a longer evolution path in In (2^, the longer parton lifetime and, as discussed in connection with Eq. (fTOl i, the enhanced 
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Figure 3: Ratios of fragmentation functions in medium and vacuum, „„j{x,Q^) /Dg yg^{x,Q^), of gluons into charged pions at = 

(100 GeV)^ where the scatterings in the medium (gluon plasma with T = 500 MeV) are included via our medium-modified evolution Eq. l|5j 
(dashed curve) and Eq. il6l (full curve). 

scattering function K for large parton energy which can give multiple scatterings that add up small energy losses to a substantial 
effect. 

The magnitude of this effect is checked in Fig. [3]using the alternative form of the scattering function in Eq. ( fTSI l. It is reassuring 
to see that the behavior as a function of x is overall the same, although somewhat different in the details. The difference between 
the two curves stems essentially from the Taylor expansion in Eq. ( fT6] l as well as the different allowed ranges for the momentum 
transfer squared \t\. As expected from this equation, the effect of the scattering becomes large at x ^ 1. 

Remembering that the quark-gluon plasma is here approximated by a stationary medium of average temperature T = 500 
MeV, one may worry that taking into account the expansion and cool-down of the plasma might lead to a smaller collisional 
energy-loss effect. This need not be the case, since the multiple scattering that can occur even in a short, early time interval when 
the plasma is denser and hotter than average, has an opposite effect. 

Jet multiplicities can be simply obtained by integrating the fragmentation function D{x,Q^) overx, but we do not consider it 
meaningful to give numerical results on this without including the suppression of soft gluon radiation. This can be done, as in 
the vacuum case, by replacing the gluon distribution D{z/u,Q^) in the DGLAP evolution term by D{z/u,u^Q^) with a reduced 
virtuality scale |28]. This will be studied in a forthcoming publication. 

V. CONCLUDING DISCUSSION 

To summarize, we have developed a formalism that combines the conventional parton cascade with parton scattering in a 
plasma of varying density and temperature. In order to give first numerical results, we simplified by assuming a plasma of con- 
stant temperature and density. This clearly demonstrated that collisional energy loss contributes substantially to jet quenching. 
The formalism developed can be included in a complete framework, where a proper expansion of the quark-gluon plasma as 
well as the momentum distributions of hard-scattered partons are included, e.g. in a Monte Carlo simulation program |22]. This 
would allow a detailed comparison to RHIC data on jet quenching. Already at this stage, it seems clear that although the overall 
behaviour is in agreement with the observed one, the magnitude is smaller than observed at RHIC. This leaves room for the 
expected additional medium-induced radiation li9i [Tol[TTIl . 

For RHIC conditions, with initial parton energy 10 GeV < £ < 20 GeV, the jet quenching effect obtained here does not 
vary strongly with this energy but depends essentially only on the momentum fraction. This gives a quenching effect which is 
essentially independent of the transverse momentum of the hadrons and thereby a nuclear modification factor R^a which is flat 
over a substantial px-interval, just as observed in the RHIC data. 

For heavy ion collisions at LHC we find a substantially increased jet quenching effect at higher jet energies due to the increased 
range of allowed momentum transfers, an increased evolution path in In Q^, the longer lifetime of the parton as well as of the 
plasma, and enhanced scattering probability. Multiple scatterings may then add up many relatively small energy losses to a 
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substantial effect. Thus, we expect, e.g., the nuclear modification factor R^a to be further reduced at the higher available at 
LHC. 

A comparison of jet profiles of heavy quarks with gluon jets may help to fix the form of the plasma scattering cross section 
more precisely. Since the momentum transfers extend to low momenta o= T, a purely perturbative calculation of the cross section 
may be questionable and an improved treatment needed 1,25,1 . Especially one has to study further how the temperature scale T 
can modify the leading behaviour. 

A recent calculation of jet evolution in the modified leading log approximation fi^ has produced differential multiplicity 
distributions. The advantage of our calculation is that it takes into account the scattering term explicitly and therefore gives 
results which depend on the plasma density and temperature. Taking transverse momentum into account explicitly as an extra 
variable, the equation could in principle be used to investigate the /Jx broadening of the parton in the medium, which is highly 
relevant for the size of the jet cone in nucleus-nucleus collisions and will be studied further 

Another issue to be studied is whether the plasma affects the non-perturbative hadronisation processes and thereby changes 
the fragmentation function. This may depend on whether partons are first forming preconfined states [S^,!!!!], as is theoretically 
supported e.g. in the large A^^ approximation, or directly form hadrons (including resonances). Intermediate prehadrons below 
the scale Qq may suffer more scatterings than final hadrons and thereby modify the hadronic spectra of heavy ion collisions. The 
non-perturbative phase after evolution down io Q = V2 GeV should be relatively more important at the lower energy scale of 
RHIC. At LHC, the larger energy implies that the formation of prehadrons and hadrons occurs later in lab frame, in particular 
for leading hadrons with a large Lorentz 7 factor We therefore expect that medium-modified hadronisation effects should be 
less and thereby our theoretical approach closer to reality. 

To conclude, we have developed a formalism to account for energy losses due to scattering in a plasma. With the indicated 
further developments this should facilitate systematic comparisons with the observed jet quenching in heavy ion collisions, 
which should contribute to a demonstration of the existence of the quark-gluon plasma and an exploration of its properties. 
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